The far infra-red SEDs of main sequence and starburst galaxies

We compare observed far infra-red/sub-millimetre (FIR/sub-mm) galaxy spectral energy distributions (SEDs) of massive galaxies ($M_{\star}\gtrsim10^{10}$ $h^{-1}$M$_{\odot}$) derived through a stacking analysis with predictions from a new model of galaxy formation. The FIR SEDs of the model galaxies are calculated using a self-consistent model for the absorption and re-emission of radiation by interstellar dust based on radiative transfer calculations and global energy balance arguments. Galaxies are selected based on their position on the specific star formation rate (sSFR) - stellar mass ($M_{\star}$) plane. We identify a main sequence of star-forming galaxies in the model, i.e. a well defined relationship between sSFR and $M_\star$, up to redshift $z\sim6$. The scatter of this relationship evolves such that it is generally larger at higher stellar masses and higher redshifts. There is remarkable agreement between the predicted and observed average SEDs across a broad range of redshifts ($0.5\lesssim z\lesssim4$) for galaxies on the main sequence. However, the agreement is less good for starburst galaxies at $z\gtrsim2$, selected here to have elevated sSFRs$>10\times$ the main sequence value. We find that the predicted average SEDs are robust to changing the parameters of our dust model within physically plausible values. We also show that the dust temperature evolution of main sequence galaxies in the model is driven by star formation on the main sequence being more burst-dominated at higher redshifts.Comment: 20 pages, 13 figures. Accepted to MNRA

G.A.S. II: Dust extinction in galaxies; Luminosity functions and InfraRed eXcess

19 pages, 18 figures, Accepted by A&AInternational audienceDust is a crucial component of the interstellar medium of galaxies. The presence of dust strongly affects the light produced by stars within a galaxy. As these photons are our main information vector to explore the stellar mass assembly and therefore understand a galaxy's evolution, modeling the luminous properties of galaxies and taking into account the impact of the dust is a fundamental challenge for semi-analytical models.We present the complete prescription of dust attenuation implemented in the new semi-analytical model: G.A.S. This model is based on a two-phase medium originating from a physically motivated turbulent model of gas structuring (G.A.S. I paper). Dust impact is treated by taking into account three dust components: Polycyclic Aromatic Hydrocarbons, Very Small Grains, and Big Grains. All three components evolve in both a diffuse and a fragmented/dense gas phase. Each phase has its own stars, dust content and geometry. Dust content evolves according to the metallicity of it associated phase.The G.A.S. model is used to predict both the UV and the IR luminosity functions from $z=9.0$ to $z=0.1$. Our two-phase ISM prescription catches very well the evolution of UV and IR luminosity functions. We note a small overproduction of the IR luminosity at low redshift ($z<0.5$). We also focus on the Infrared-Excess (IRX) and explore its dependency with the stellar mass, UV slope, stellar age, metallicity and slope of the attenuation curves. Our model predicts large scatters for relations based on IRX, especially for the IRX-$\beta$ relation. Our analysis reveals that the slope of the attenuation curve is more driven by absolute attenuation in the FUV band than by disk inclination. We confirm that the age of the stellar population and the slope of the attenuation curve can both shift galaxies below the fiducial star-birth relation in the IRX-$\beta$ diagram

Star formation properties of sub-mJy radio sources

We investigate the star formation properties of ~800 sources detected in one of the deepest radio surveys at 1.4 GHz. Our sample spans a wide redshift range (~0.1 - 4) and about four orders of magnitude in star formation rate (SFR). It includes both star forming galaxies (SFGs) and active galactic nuclei (AGNs), further divided into radio-quiet and radio-loud objects. We compare the SFR derived from the far infrared luminosity, as traced by Herschel, with the SFR computed from their radio emission. We find that the radio power is a good SFR tracer not only for pure SFGs but also in the host galaxies of RQ AGNs, with no significant deviation with redshift or specific SFR. Moreover, we quantify the contribution of the starburst activity in the SFGs population and the occurrence of AGNs in sources with different level of star formation. Finally we discuss the possibility of using deep radio survey as a tool to study the cosmic star formation history.Comment: 18 pages, 14 figures, 1 table (available in its entirety as ancillary data

ALMA observations of atomic carbon in z~4 dusty star-forming galaxies

We present ALMA [CI]($1-0$) (rest frequency 492 GHz) observations for a sample of 13 strongly-lensed dusty star-forming galaxies originally discovered at 1.4mm in a blank-field survey by the South Pole Telescope. We compare these new data with available [CI] observations from the literature, allowing a study of the ISM properties of $\sim 30$ extreme dusty star-forming galaxies spanning a redshift range $2 < z < 5$. Using the [CI] line as a tracer of the molecular ISM, we find a mean molecular gas mass for SPT-DSFGs of $6.6 \times 10^{10}$ M$_{\odot}$. This is in tension with gas masses derived via low-$J$ $^{12}$CO and dust masses; bringing the estimates into accordance requires either (a) an elevated CO-to-H$_2$ conversion factor for our sample of $\alpha_{\rm CO} \sim 2.5$ and a gas-to-dust ratio $\sim200$, or (b) an high carbon abundance $X_{\rm CI} \sim 7\times10^{-5}$. Using observations of a range of additional atomic and molecular lines (including [CI], [CII], and multiple transitions of CO), we use a modern Photodissociation Region code (3D-PDR) to assess the physical conditions (including the density, UV radiation field strength, and gas temperature) within the ISM of the DSFGs in our sample. We find that the ISM within our DSFGs is characterised by dense gas permeated by strong UV fields. We note that previous efforts to characterise PDR regions in DSFGs may have significantly underestimated the density of the ISM. Combined, our analysis suggests that the ISM of extreme dusty starbursts at high redshift consists of dense, carbon-rich gas not directly comparable to the ISM of starbursts in the local Universe.Comment: 21 pages, 12 figures. Accepted for publication in MNRA

Sub-kiloparsec Imaging of Cool Molecular Gas in Two Strongly Lensed Dusty, Star-Forming Galaxies

We present spatially-resolved imaging obtained with the Australia Telescope Compact Array (ATCA) of three CO lines in two high-redshift gravitationally lensed dusty star-forming galaxies, discovered by the South Pole Telescope. Strong lensing allows us to probe the structure and dynamics of the molecular gas in these two objects, at z=2.78 and z=5.66, with effective source-plane resolution of less than 1kpc. We model the lensed emission from multiple CO transitions and the dust continuum in a consistent manner, finding that the cold molecular gas as traced by low-J CO always has a larger half-light radius than the 870um dust continuum emission. This size difference leads to up to 50% differences in the magnification factor for the cold gas compared to dust. In the z=2.78 galaxy, these CO observations confirm that the background source is undergoing a major merger, while the velocity field of the other source is more complex. We use the ATCA CO observations and comparable resolution Atacama Large Millimeter/submillimeter Array dust continuum imaging of the same objects to constrain the CO-H_2 conversion factor with three different procedures, finding good agreement between the methods and values consistent with those found for rapidly star-forming systems. We discuss these galaxies in the context of the star formation - gas mass surface density relation, noting that the change in emitting area with observed CO transition must be accounted for when comparing high-redshift galaxies to their lower redshift counterparts.Comment: 14 pages, 7 figures; accepted for publication in Ap

SPT0346-52: Negligible AGN Activity in a Compact, Hyper-starburst Galaxy at z = 5.7

We present Chandra ACIS-S and ATCA radio continuum observations of the strongly lensed dusty, star-forming galaxy SPT-S J034640-5204.9 (hereafter SPT0346-52) at $z$ = 5.656. This galaxy has also been observed with ALMA, HST, Spitzer, Herschel, APEX, and the VLT. Previous observations indicate that if the infrared (IR) emission is driven by star formation, then the inferred lensing-corrected star formation rate ($\sim$ 4500 M_{\sun} yr$^{-1}$) and star formation rate surface density $\Sigma_{\rm SFR}$ ($\sim$ 2000 M_{\sun} {yr^{-1}} {kpc^{-2}}) are both exceptionally high. It remained unclear from the previous data, however, whether a central active galactic nucleus (AGN) contributes appreciably to the IR luminosity. The {\it Chandra} upper limit shows that SPT0346-52 is consistent with being star-formation dominated in the X-ray, and any AGN contribution to the IR emission is negligible. The ATCA radio continuum upper limits are also consistent with the FIR-to-radio correlation for star-forming galaxies with no indication of an additional AGN contribution. The observed prodigious intrinsic IR luminosity of (3.6 $\pm$ 0.3) $\times$ 10$^{13}$ L_{\sun} originates almost solely from vigorous star formation activity. With an intrinsic source size of 0.61 $\pm$ 0.03 kpc, SPT0346-52 is confirmed to have one of the highest $\Sigma_{SFR}$ of any known galaxy. This high $\Sigma_{SFR}$, which approaches the Eddington limit for a radiation pressure supported starburst, may be explained by a combination of very high star formation efficiency and gas fraction.Comment: 8 pages, 6 figures, accepted for publication in Ap

THE EVOLVING INTERSTELLAR MEDIUM OF STAR-FORMING GALAXIES SINCE z=2 AS PROBED BY THEIR INFRARED SPECTRAL ENERGY DISTRIBUTIONS

Using data from the mid-infrared to millimeter wavelengths for individual galaxies and for stacked ensembles at 0.5 1012 L ☉). For galaxies within the MS, we show that the variations of specific star formation rates (sSFRs = SFR/M *) are driven by varying gas fractions. For relatively massive galaxies like those in our samples, we show that the hardness of the radiation field, langUrang, which is proportional to the dust-mass-weighted luminosity (L IR/M dust) and the primary parameter defining the shape of the IR spectral energy distribution (SED), is equivalent to SFE/Z. For MS galaxies with stellar mass log (M */M ☉) ≥ 9.7 we measure this quantity, langUrang, showing that it does not depend significantly on either the stellar mass or the sSFR. This is explained as a simple consequence of the existing correlations between SFR-M *, M *-Z, and M gas-SFR. Instead, we show that langUrang (or equally L IR/M dust) does evolve, with MS galaxies having harder radiation fields and thus warmer temperatures as redshift increases from z = 0 to 2, a trend that can also be understood based on the redshift evolution of the M *-Z and SFR-M * relations. These results motivate the construction of a universal set of SED templates for MS galaxies that are independent of their sSFR or M * but vary as a function of redshift with only one parameter, langUrang

The Cosmic Far-Infrared Background Buildup Since Redshift 2 at 70 and 160 microns in the COSMOS and GOODS fields

The Cosmic Far-Infrared Background (CIB) at wavelengths around 160 {\mu}m corresponds to the peak intensity of the whole Extragalactic Background Light, which is being measured with increasing accuracy. However, the build up of the CIB emission as a function of redshift, is still not well known. Our goal is to measure the CIB history at 70 {\mu}m and 160 {\mu}m at different redshifts, and provide constraints for infrared galaxy evolution models. We use complete deep Spitzer 24 {\mu}m catalogs down to about 80 {\mu}Jy, with spectroscopic and photometric redshifts identifications, from the GOODS and COSMOS deep infrared surveys covering 2 square degrees total. After cleaning the Spitzer/MIPS 70 {\mu}m and 160 {\mu}m maps from detected sources, we stacked the far-IR images at the positions of the 24 {\mu}m sources in different redshift bins. We measured the contribution of each stacked source to the total 70 and 160 {\mu}m light, and compare with model predictions and recent far-IR measurements made with Herschel/PACS on smaller fields. We have detected components of the 70 and 160 {\mu}m backgrounds in different redshift bins up to z ~ 2. The contribution to the CIB is maximum at 0.3 <= z <= 0.9 at 160{\mu}m (and z <= 0.5 at 70 {\mu}m). A total of 81% (74%) of the 70 (160) {\mu}m background was emitted at z < 1. We estimate that the AGN relative contribution to the far-IR CIB is less than about 10% at z < 1.5. We provide a comprehensive view of the CIB buildup at 24, 70, 100, 160 {\mu}m. IR galaxy models predicting a major contribution to the CIB at z < 1 are in agreement with our measurements, while our results discard other models that predict a peak of the background at higher redshifts. Our results are available online http://www.ias.u-psud.fr/irgalaxies/ .Comment: Accepted in Astronomy & Astrophysic

The Rest-Frame Submillimeter Spectrum of High-Redshift, Dusty, Star-Forming Galaxies

We present the average rest-frame spectrum of high-redshift dusty, star-forming galaxies from 250-770GHz. This spectrum was constructed by stacking ALMA 3mm spectra of 22 such sources discovered by the South Pole Telescope and spanning z=2.0-5.7. In addition to multiple bright spectral features of 12CO, [CI], and H2O, we also detect several faint transitions of 13CO, HCN, HNC, HCO+, and CN, and use the observed line strengths to characterize the typical properties of the interstellar medium of these high-redshift starburst galaxies. We find that the 13CO brightness in these objects is comparable to that of the only other z>2 star-forming galaxy in which 13CO has been observed. We show that the emission from the high-critical density molecules HCN, HNC, HCO+, and CN is consistent with a warm, dense medium with T_kin ~ 55K and n_H2 >~ 10^5.5 cm^-3. High molecular hydrogen densities are required to reproduce the observed line ratios, and we demonstrate that alternatives to purely collisional excitation are unlikely to be significant for the bulk of these systems. We quantify the average emission from several species with no individually detected transitions, and find emission from the hydride CH and the linear molecule CCH for the first time at high redshift, indicating that these molecules may be powerful probes of interstellar chemistry in high-redshift systems. These observations represent the first constraints on many molecular species with rest-frame transitions from 0.4-1.2mm in star-forming systems at high redshift, and will be invaluable in making effective use of ALMA in full science operations.Comment: 19 pages, 10 figures (2 in appendices); accepted for publication in Ap

Merger driven star-formation activity in Cl J1449+0856 at z=1.99 as seen by ALMA and JVLA

We use ALMA and JVLA observations of the galaxy cluster Cl J1449+0856 at z=1.99, in order to study how dust-obscured star-formation, ISM content and AGN activity are linked to environment and galaxy interactions during the crucial phase of high-z cluster assembly. We present detections of multiple transitions of 12CO, as well as dust continuum emission detections from 11 galaxies in the core of Cl J1449+0856. We measure the gas excitation properties, star-formation rates, gas consumption timescales and gas-to-stellar mass ratios for the galaxies. We find evidence for a large fraction of galaxies with highly-excited molecular gas, contributing >50% to the total SFR in the cluster core. We compare these results with expectations for field galaxies, and conclude that environmental influences have strongly enhanced the fraction of excited galaxies in this cluster. We find a dearth of molecular gas in the galaxies' gas reservoirs, implying a high star-formation efficiency (SFE) in the cluster core, and find short gas depletion timescales τ<0.1-0.4 Gyrs for all galaxies. Interestingly, we do not see evidence for increased specific star-formation rates (sSFRs) in the cluster galaxies, despite their high SFEs and gas excitations. We find evidence for a large number of mergers in the cluster core, contributing a large fraction of the core's total star-formation compared with expectations in the field. We conclude that the environmental impact on the galaxy excitations is linked to the high rate of galaxy mergers, interactions and active galactic nuclei in the cluster core